Electrical transport through single-molecule junctions: from molecular orbitals to 

conduction channels 
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We present an atomistic theory of electronic transport through single organic molecules that 
reproduces the important features of the current-voltage (I-V) characteristics observed in recent 
experiments. We trace these features to their origin in the electronic structure of the molecules and 
their local atomic environment. We demonstrate how conduction channels arise from the molecular 
orbitals and elucidate which specific properties of the individual orbitals determine their contribution 
to the current. 

PACS numbers: 73.40.Jn, 73.40.Cg, 73.40.Gk, 85.65.+h 
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Introduction. — Present trends in the miniaturiza- 
tion of electronic devices suggest that ultimately sin- 
gle molecules may be used as electronically active ele- 
ments in a variety of applications [ffl, ^ . Recent advances 
in the manipulation of single molecules now permit to 
contact an individual molecule between two electrodes 
(see Fig. 1) and measure its electronic transport proper- 
ties H fl Hj H @, U ■ In contrast to single-electron transis- 
tors based on metallic islands B , molecular devices have 
a more complicated, but in principle tunable, electronic 
structure. Interesting and novel effects, such as negative 
differential conductance []l0|| , were observed in some of 
these experiments, which still, by-and-large, beg theo- 
retical explanation. In addition to generic principles of 
nanoscale physics, e.g. Coulomb blockade |ll[ [l2|, the 
chemistry and geometry of the molecular junction emerge 
as the fundamental tunable characteristics of molecular 
junctions § § g £| 

In this paper we present an atomistic theory that 
bridges traditional concepts of mesoscopic and molecu- 
lar physics to describe transport through single organic 
molecules in qualitative agreement with recent break- 
junction experiments [0. We combine ab initio quantum 
chemistry calculations with non-equilibrium Green func- 
tions techniques to illustrate the emergence of conduction 
channels in a single-molecule junction from the molecular 
orbitals (MO). We further show how the specific proper- 
ties of individual MOs are reflected in their contribution 
to the current. Using this data we provide insight into 
the microscopic origin of the nonlinear I-V characteris- 
tics observed experimentally and correlate their features 
to the specific properties of the molecule. 

Our approach naturally accounts for the experimen- 
tal observations and indicates that the current in these 
molecular junctions is mainly controlled by the electronic 
structure of the molecules and their local environment. 
We demonstrate that many molecular orbitals partici- 
pate in a single conduction channel and provide exam- 
ples where, surprisingly, the current is not dominated by 
the contribution from the energetically closest MOs. The 



theory provides a quantitative criterion to judge the im- 
portance of individual MOs to the current and thus paves 
the way for the a priori design of molecular transport 
properties. 

Theoretical model. — We calculate the current through 
a single molecule attached to metallic electrodes by a gen- 
eralization of an earlier analysis of transport in atomic- 
size contacts [Q, similar in spirit to Refs. jlj, [l7|. Since 
the conductance is mainly determined by the narrowest 
part of the junction, only the electronic structure of this 
"central cluster" must be resolved in detail. It is there- 
fore sensible to decompose the overall Hamiltonian of the 
molecular junction as 



H = H L + H R + H c + V, 



(1) 



where He describes the "central cluster" of the system, 
Hl,r describe the left and right electrode respectively, 
and V gives the coupling between the electrodes and the 
central cluster (see Fig. 1). 

In this study the electronic structure of the "cen- 
tral cluster" is calculated within the density functional 
(DFT) approximation The left and right reservoirs 
are modeled as two perfect semi-infinite crystals of the 
corresponding metal using a tight-binding parameteri- 
zation |^0| . Finally, V describes the coupling between 
the leads and the central cluster and takes the form: 
^ = Ylij v ij(d-lcj + h.c.). The hopping elements u,j be- 
tween the lead orbitals d\ and the MOs of the central 
cluster 6j are obtained by reexpressing the MOs of the 
central cluster via a Lowdin transformation in terms of 
atom-like orbitals, and then using the mentioned tight- 
binding parameterization pp[ . 

The "central cluster" is not necessarily confined to the 
molecule, but may, in principle, contain arbitrary parts of 
the metallic electrode. The inclusion of part of the leads 
in the ab initio calculation was shown to improve the 
description of the molecule-leads coupling , in partic- 
ular regarding charge transfer between the molecule and 
the electrodes. The Fermi energy of the overall system 
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H C N O S Au leads 



FIG. 1: Scheme of the single-molecule contacts analyzed in 
this work. The two organic molecules attached to gold elec- 
trodes, which were experimentally investigated in Ref. || , are 
referred to as: (a) "symmetric molecule" , and (b) "asymmet- 
ric molecule". 



is determined by the charge neutrality condition of the 
central cluster. 

In order to obtain the current for a constant bias 
voltage, V, between the leads, we make use of non- 
equilibrium Green function techniques. Since the Hamil- 
tonian of Eq. (1) does not contain inelastic interactions, 
the current follows from the Landauer formula [Ell 



2<= f°° 
I=-^J_ deTr{W}[f(e 



eV/2)-f(e + eV/2)], (2) 



In order to provide a deep understanding of the electronic 
transport, we analyze the current in terms of conduction 
channels, defined as eigenfunctions of ffi . Such analy- 
sis allows to quantify the contribution to the transport 
of every individual molecular level. In our approach the 
channels arise as a linear combination of the molecular 
orbitals \<frj) of the central cluster, i.e. |c) — J2j a cj\ ( t ) j) > 
and the corresponding eigenvalues determine their contri- 
bution to the conductance. Ultimately, this information 
concerning the channels could eventually be measured 
using superconducting electrodes p3|. 




where / is the Fermi function and t is the energy and 
voltage dependent transmission matrix given by 



FIG. 2: Total density of states (TDOS) of the molecule and 
zero-bias total transmission as a function of the energy for 
both molecules. The Fermi energy is set to zero. 



i(e,V) = 2f 1 I ! 2 (e 



eV/2)G r c (e,V)f R /2 (e + eV/2). (3) 



The scattering rate matrices are given by Tl^r = 
Im(£L,i?), where are the self-energies which contain 
the information of the electronic structure of the leads 
and their coupling to the central cluster. They can be 
expressed as Tj L ..r{^) = vcl,r9l,r(£)vl..r.c, V being the 
hopping matrix which describes the connection between 
the central cluster and the leads. Ql,r are the Green 
functions of the uncoupled leads (semi- infinite crystals), 
which are calculated using decimation techniques P2| . 
The Green functions of the central cluster are given by 



G c (e, V)= el -He- t L {e - eV/2) -t R (e + eV/2) 



(4) 



Results and discussions. — We now use the method de- 
scribed above to analyze the experiment of Reichert et 
al. H . The two organic molecules investigated are shown 
in Fig. 1, and as indicated in the caption, they will be re- 
ferred to as symmetric and asymmetric molecule. For the 
description of the gold reservoirs, we use a basis with the 
atomic-like 5d, 6d, 6p orbitals, and for the central clus- 
ter we use the LANL2DZ basis Q for all atoms. The 
DFT calculations were performed with the Becke three- 
parameter hybrid functional using the Lee, Yang and 
Parr correlation functional ]25| at zero field (^6|. In the 
calculations reported here one additional gold atom was 
included on either side of the molecule. Experimentally, 
both molecules were contacted several times and the na- 
ture of the I-V characteristics was found to vary with 
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FIG. 3: (a-d) Charge-density plots of four molecular orbitals of the central cluster for the symmetric molecule. Panel (a) 
displays the HOMO and (c) the LUMO, which is twofold degenerate, (b) shows a confined orbital that contributes little to the 
current, while the MO in (d) is almost as important as the LUMO despite its difference in energy. Panel (e) shows the total 
density of states of the central cluster (dotted line) and the individual contributions of the four molecular orbitals (color lines). 
The level positions are indicated on top of this panel. The contributions of the different MOs to the conduction channel at the 
Fermi energy (set to zero) are: \a a \ 2 = 0.007, |a b | 2 = 10" 11 , ja c | 2 = 0.06, ja d | 2 = 0.02. 



the quality of the contact. For this reason, theory can 
presently aim to elucidate important reproducible fea- 
tures of the experiment under the assumption that the 
contact to the electrodes is well defined. Since there is 
no direct experimental information regarding the geom- 
etry of the molecule and its attachment to the leads, the 
overall geometry of the central cluster was relaxed with- 
out additional constraints in our calculations, resulting 
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FIG. 4: I-V characteristics and differential conductance for 
the symmetric (a) and asymmetric (b) molecules. The inset 
in (b) shows the charge-density plot of the HOMO for the 
asymmetric molecule. Notice the intrinsic asymmetry of the 
charge distribution in the gold atoms. 



in the Au atom being out of the molecular plane. 

Let us start by analyzing the linear response regime. 
In Fig. 2 we show for both cases the total density of states 
(TDOS) of the molecule and the zero-bias total transmis- 
sion as a function of energy. As can be seen in the TDOS, 
in both cases the covalent bond between Au and S results 
in a strong hybridization between the molecular orbitals 
and the extended states of the metallic electrodes. The 
formation of wide energy bands and the disappearance of 
the gaps of the discrete molecular spectrum suggest the 
absence of Coulomb blockade in this type of molecular 
junctions. 

The zero-bias total transmission as a function of en- 
ergy follows closely the TDOS. The transmission is dom- 
inated overwhelmingly by a single channel in the energy 
window shown in Fig. 2, and the corresponding eigen- 
values of ffi at the Fermi energy |27j are T sym = 0.014 
and T asyrn = 0.006. The decomposition of this channel 
into molecular orbitals provides us information on the 
relevance of the different molecular levels. Fig. 3 (a)-(d) 
show charge-density plots for some representative MOs 
and (e) shows their individual contribution to the TDOS. 
Fig. 3a shows that the highest occupied molecular or- 
bital (HOMO) is confined to the interior of the molecule 
and its weight at the gold atoms is rather small. Con- 
sequently, in spite of its privileged energy position, the 
HOMO does not give a significant contribution to the cur- 
rent. The lowest unoccupied molecular orbital (LUMO), 
see Fig. 3c, exhibits the opposite behavior, i.e. it is very 
well coupled to the leads through the 6s atomic orbital 
of the gold atoms (notice that it has width of about 4 eV 
in the density of states), but the charge is mainly local- 
ized on the Au and S atoms. The interplay of these two 
factors yields a contribution of w 6% of the total current. 
Fig. 3(b,d) shows two further MOs with similar energy 
but very different contribution to the channel. While the 
localized MO (b) carries almost no current, the extended 
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and well coupled MO (d) has significant weight. Conse- 
quently there are three ingredients which determine the 
contribution of a MO to the current: (i) its energy po- 
sition (distance to the Fermi energy), (ii) its bridging 
extent (whether it is extended or localized), and (iii) its 
coupling to the leads. Our analysis provides a counterex- 
ample to the conventional wisdom that the HOMO and 
the LUMO dominate the transport properties. 

Fig. 4 shows the I-V curves for both molecules in the 
voltage range investigated in Ref. B. Both the order 
of magnitude and shape of the current and conductance 
agree qualitatively with the experimental results. There 
is no pronounced voltage dependence of the transmis- 
sion due to the smooth density of states of the gold 
electrodes in the energy region explored here. The non- 
linearities in these I-V curves can be then understood by 
a simple inspection of the energy dependence of the zero- 
bias transmission. For instance, the pronounced increase 
in the conductance around 1 V is due to the fact that 
we approach the resonant condition for the HOMO and 
LUMO. 

In agreement with the experiment, the I-V of the sym- 
metric molecule is symmetric with respect to voltage 
inversion, while the one of the asymmetric molecule is 
asymmetric. According to Eq. 2 we note that the asym- 
metry of the MO is not a sufficient condition for an asym- 
metric I-V characteristic, as evidenced by the experimen- 
tal I-V near zero bias. Only the energy dependence of the 
left/right scattering rates in Eq. 3 induces an asymme- 
try into the I-V characteristic. Such asymmetries can 
arise either from differences in the local atomic struc- 
ture of the electrodes or from intrinsic properties of the 
molecule. The latter arise from the asymmetric charge 
distribution of the MOs (see inset Fig. 4b), resulting in 
different couplings to the leads. Even for the symmet- 
ric molecule, we were able to induce asymmetries into 
the I-V characteristic by distorting the geometry of one 
of the lead fragments in the central cluster. This fact 
was nicely demonstrated in the experiment (see Fig. 5 
in Ref. (8|). We investigated several other scenarios re- 
garding the number of gold atoms, their geometry and 
the coupling and found predictable variations of the I- 
V's with theses changes. If, for example, the gold-sulfur 
bond is stretched by about 5%, the current changes by 
about 70%, but the position of the peaks in the conduc- 
tance and their peak/valley ratios are only marginally 
affected. 

Conclusions. — We have presented an atomistic semi- 
quantitative description of non-linear transport through 
a single molecule junction. We were able to attribute dis- 
tinctive features of the I-V's of the symmetric and asym- 
metric molecule to their individual molecular levels ob- 



tained from ab initio calculations. The resolution of con- 
ductance into conduction channels permits an analysis of 
the contributions of individual orbitals to overall trans- 
port. In this way we provide an understanding which 
can be valuable for the future engineering of molecular 
devices. 
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